Charge Qubit Purification by an Electronic Feedback Loop 
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We propose the manipulation of an isolated qubit by a simple instantaneous closed-loop feedback 
scheme in which a time-dependent electronic detector current is directly back-coupled into qubit 
parameters. As specific detector model we employ a capacitively coupled single-electron transistor. 
We demonstrate the stabilization of pure delocalized qubit states above a critical detector-qubit 
coupling. This electronic purification is independent of the initial qubit state and is accomplished 
after few electron jumps through the detector. Our simple scheme can be used for the efficient 
and robust initialization of solid-state qubits in quantum computational algorithms at arbitrary 
temperatures. 

PACS numbers: 03.65.Ta, 05.40.-a, 73.23.-b, 85.35.-p 



The ability to measure and control individual quantum 
systems (e.g. qubits) is a major requirement for contem- 
porary and future information technologies [l| . The mea- 
surement process - the contact between classical mea- 
surement device and nano-scale object - plays a key role 
in controlling quantum systems and mostly introduces 
a source of dccohcrcncc. State-of-the-art measurement 
technologies in solid-state environments mainly rely on 
the change of current flow through an on-chip detector 
due to Coulomb repulsion with the charges in the mea- 
sured system (see e.g. Refs. The continuous mon- 
itoring of charge qubits by such charge detectors yields 
relaxation and decoherence as studied experimentally e.g. 
in Rcf. Q and theoretically in Refs. However, the de- 
tector signal contains detailed information on the qubit 
dynamics which can be used for its control. 

An elegant control strategy, which has been success- 
fully applied in a wide variety of physical and engineering 
systems Q , is stabilization of favored system states due 
to feedback. For quantum systems such as qubits, sev- 
eral theoretical proposals to enhance quantum coherence 
by feedback have been made. These are mainly based 
upon sophisticated techniques, e.g., the on-chip conver- 
sion of detector outputs into qubit density matrix evolu- 
tions 0, 01 or monitoring the phase of qubit oscillations 
by detector current quadratures However, a much 
simpler scheme was recently used in the first experimen- 
tal use of feedback to increase quantum coherence of a 
spin qubit 

Here, we study a charge qubit consisting of a single 
electron confined in a two-state system which is continu- 
ously monitored by the current through a single-electron 
transistor (SET) [see Fig. [1]. The SET current wiU be 
directly fed back into qubit parameters which allow the 
stabilization of arbitrary pure qubit states. 
Model. - The SET electronic level is shifted depending 
on where the electron is located in the two-level system 
due to Coulomb interaction. Moreover, the charge qubit 
electron is assumed to be coupled to thermal phonons of 
the surrounding matrix material. The full Hamiltonian 



reads 
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HsET + H^ + -dU[Ua, + {Ut + Ub)t], (1) 
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with e = Et — Sb and Aq cl-q -r uq. 
with Hamiltonian Hset = ^dd'' d + J2ka=L,Rl'^kacl^Cka + 
{tkaCkad^ + h.c.)] coutaius a single electronic level 
(with creation/annihilation operators d^/d) which is 
electrostatically coupled to the charge qubit levels by 
U = Ut — Ub > 0- The charge qubit Hamiltonian is 
given in the Bloch representation with Pauli matrices 
= \t) (t\ - \b) {b\ and (7^ = \t) {b\ + \b) {t\. Thephononic 
modes with wave vector Q and energy Qq are diagonally 
coupled to the qubit assuming weak coupling Tc- 
Equations of Motion. - The standard Born-Markov- 
secular approximation yields a quantum master equa- 
tion p{t) = Cp{t) for the system (SET -I- qubit) den- 
sity matrix p{t) with the Liouvillian superoperator C 
where the SET leads and the phonon environment are 
treated in thermal equilibrium. Due to the secular ap- 
proximation the 4 populations and 2 coherences con- 
tained in p{t) decouple in the energy eigenbasis. In 
the Markovian limit, tunneling rates and bath expecta- 
tion values are all evaluated at the transition frequen- 
cies of the system. Since the SET Coulomb shift de- 
pends on the position of the charge qubit electron, this 
introduces for finite SET bias effective tunneling rates 
Ta and Fjj for lower and upper electron position, re- 
spectively, such that the (time and ensemble averaged) 
SET current comprises damped oscillations for FJ^ 7^ Fq, 
see Fig. [Hd) . These are caused by the qubit evolution 
[Fig. Hie)]. With the SET occupation n{t) = Tr[np(<)] 
and the Bloch vector components conditioned on whether 
the SET is occupied, (t) = Tr[(n (g) CTa)p(0] and 
S'q (t) = Tr{[(l — fi) (g) d-a]p{t)} {a = x, y, z) respectively, 
such that S°'{t) = Sg{t) + S^{t), the quantum master 
equation has the following convenient vector form (omit- 
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FIG. 1: (color online), (a) Sketch of the device: charge qubit 
and single-electron transistor (SET) as capacitively coupled 
detector; a quantum point contact (QPC) is attached to the 
SET to resolve lL{t) and lR{t) for the closed-loop feedback, 
(b) Single realization of the SET occupation n{t) measured 
by the QPC current /Qpc(i)- (c) Single realizations of SET 
currents /^(i) (blue, solid) and lR{t) (dashed, green), (d) 
SET current {I){t), ensemble averaged (.) and time averaged 
/ = \ ^l^^ I{t')dt' with I(t) = \[lL{t) + lR{t)], exhibits 
damped oscillations caused by the coupling to the qubit. (e) 
Bloch vector component ^^(t) (qubit occupation): ensemble 
average (red, dashed) shows oscillations damped by the pres- 
ence of the SET, single realization (blue, solid) exhibits free 
evolution interrupted whenever an electron jumps at the SET. 
Parameters: F'/F = 4, Tc/F = 1, e = 0, f//F = 0.5, measure- 
ment time Fr = 0.01, initial state: S"'(0) = 1, n(0) = 0. 



ting the time dependence): 

50 = Ao X So+I?LSo + Si?Si+non+(77^,0,7^)'^, 

51 = Ai X Si +PhSi+6lSo + nm- (0,0, 7^)^,(2) 

with s,(t) ^ (5f (t), sy{t),st{t)Y, 7^ ^Wa± r;), 

7± = 7^+71, no = -ni = -(77^,0,7^)'^, and 

"7a+ - Vp 
-7„+ 

Ba = Diag (v/TJ^, yrj^, 7+) , (3) 



with the phonon rates r]p and r]r = 77+ — explicitly 
defined in Ref. [l^l- For zero qubit bias (e = 0) they 
are rjr = gTTTce~^'^''^^'= and rjp = 77^ coth (/JTc) with the 
dimensionless coupling strength g, the cut-off frequency 
LUc, and the inverse temperature /3. The free dynamics 
of Sj{t) in Eqs. ^ can be considered as a precession 

around vector Aj = —(2Tc,0,e + jU)'^ with frequency 
A = y/e"^ + 4r2 (j = 0) or A' = v^(e + Uy + (j = 
1), respectively. The effect of dissipation is contained 
in Va'- damping of all Bloch vector components by the 
presence of the SET with 7^ or 7]^ and by the electron- 
phonon coupling which leads to an additional decay and 
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FIG. 2: (color online). Bloch vector component S~{t) vs. 
time rt and feedback parameter ©ij^ (a) and Qrt (b); (c) 
S'^{t) cross-sections at Ors — —1 (blue, dash-dotted), Qrt = 
—0.5 (red, full), 0.5 (green, dashed). Parameters: e = 0, 
U/r = 0.5, F'/F = 3, Tc/F = 1. Initial state is the complete 
statistical mixture with empty SET. With feedback damped 
qubit oscillations can be induced which would not occur for 
vanishing SET-qubit coupling F = F'. 



scattering of the S'f^^(i). The conditioned Bloch vec- 
tors Sj {t) are mutually coupled via the diagonal matrices 
B]^/ji, i.e., by electron hopping onto the SET and out of 
it, respectively. We provide an analytical solution in the 
online supplement The SET occupation couples to 
the qubit and is not affected by its dynamics for Fq = Fjj . 
We note that the Coulomb interaction U enters Eqs. ([2]) 
implicitly in F' and explicitly in Ai. Whenever we set 
U = m the following we do not exclude T' y^T in order 
to maintain the coupling but rather neglect the physical 
back-action on the qubit. 

Detector counting statistics and active measurement. - 
The information about the number of charges Ua trans- 
ferred through a multi-lead conductor (our SET detec- 
tor) at the a-th lead can be gained by the help of tt.^- 
rcsolved quantum master equations (see e.g. Ref. [l^ for 
detailed discussions). The unraveled Fourier-transformed 
master equation for unidirectional transport is p = [Cq + 
Y.a (^'^'"J'a]p = C{{xa})p, whcrc the Ja{£o) denote jump 
(non-jump) operators and Xa are the counting fields and 
the system state p is not actively controlled (passive 
measurement). In contrast, when active control oper- 
ations on the system depending on the measurement 
outcomes arc performed after single jump events, the 
Fourier-transformed master equation generalizes to 



^p({xJ,i) = 



Co 



p{{Xa},t) (4) 



where the ]Ca describe the control operation In con- 
trast to passive measurements, in an active (control) 
measurement the (scalar) counting fields Xa are thus up- 
graded to superoperators, ixa — > iXa + fCa- 
Feedback scheme. - We choose the feedback of the SET 
current acting on Tc or/and on e which are accessi- 
ble in present-days setups by appropriate gate voltages. 
The corresponding feedback Hamiltonian is defined by 
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^fb = J2a=L/R^a{t)ha with ha = Qae^z + QaT^x and 

the feedback parameters OqT and Qae- Experimentally, 
single jump events I]^/fi{t) may neither be resolved nor 
measured independently in the SET circuit. Therefore, 
we propose an additional quantum-point contact (QPC) 
weakly attached to the SET as sketched in Fig. HJa). 
Its telegraph-like current signal Iqpcit) provides the oc- 
cupation of the SET [Fig. Wlh)]. For F'^ > Fa the 
upwards/downwards steps are related to single jumps 
into/out of the SET and, consequently, lL{t)/lR{t) are 
trains of (5-kicks [Fig. [D^c)]. The supcroperators /C^/j^ 
in the feedback master equation (|4|) are then given by 
K-aP = —i[ha,p]- With feedback into s/Tc, the diagonal 
matrices B]^ //j 131 become rotation matrices around the 
qubit x/z-axis so that whenever an electron jumps 
into/out of the SET the qubit state is rotated by the 
corresponding angle. 

Feedback-induced qubit oscillations. - As a remarkable 
effect of this Markovian feedback, we find damped qubit 
oscillations even when the qubit is initially prepared in 
a completely mixed state (i.e., in the origin of the Bloch 
sphere: S = 0). For low clcctron-phonon coupling with- 
out feedback or for vanishing SET-qubit coupling 
(F' = F, i.e., the qubit kicked by a random sequence 
of (5-pulses) the qubit would remain completely mixed. 
Fig. [2] shows the oscillations in (t) for a feedback from 
IfiX^) [^L{t) provides similar results] into (a) e and (b) 
Tc- The feedback into £ [Fig. EJa)] induces oscillations 
with the frequency ujq = \y/{8Tc)'^ — 7I (e = 0) for 
8Tc > 7s where 75 ^{Vt — VT')'^ is the damping 
due to SET-qubit coupling [111 and additional feedback 
damping which increases with Qrc- At large times, the 
qubit state is again completely mixed. In contrast, for 
the Tc— feedback [Fig. ^b)] a frequency shift towards 
highcr(lower) frequencies for Qrt > (<) occurs, and in 
the long-term limit a localization [S^{t — > 00) 7^ 0] takes 
place. The transient oscillations of the qubit occupation 
[Fig. [2|c)] can be directly measured as damped current 
oscillations [compare Figs.[T{d) and (e)]. 

Noise spectrum. - The total noise spectrum of the 
SET current is computed by (Ramo-Shockley theorem) 
S{uj) ^ jJ2a b=L B. ^ab{^) for a Symmetric SET with 
Sabiuj) = Tr{ [Xfc + Ja^Wb + Jani-uj)Jb]p} whcrc 
n{uj) = [-iwl-/:(0,0)]-\ steady state p with £(0,0)p = 
0, and the jump operators Ja = di^^C{xL,XR)\(Q,o)^ 
Ja,b = di^a9ixi^XL,XR)\(o,o)- The spectrum provides 
detailed information on the qubit dynamics such as os- 
cillation frequencies and relaxation rates. Without feed- 
back, S(lli) shows peaks at the bare qubit oscillation fre- 
quency A (for 7s <C Tc) and at the shifted frequency due 
to Coulomb interaction A' for nonvanishing and weak 
SET-qubit coupling F'/F ^ 1, see inset of Fig. EKa). We 
can further use it to monitor the effect of feedback on 
the qubit dynamics, as shown in Fig. |3Ka). The most in- 
teresting feedback scenario is when both /l(0 and /i?(t) 
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FIG. 3: (color online), (a) SET noise spectrum S(cj) vs. uj 
for various feedback parameter Ot = Qlt = Ort = —0.4 
(dashed, red), —0.2 (dotted, green), (thick solid, black), 
0.2 (dash-dotted, blue), 0.4 (thin solid, violet); (Inset) S{uj) 
for Qt = 0, r'/r = 1.5, U/r = 5. (b) single realization of 
Bloch vector components S^''^{t) and SET occupation n{t) vs. 
time Ft with initially mixed state S = 0; ensemble averaged 
quantities (dashed lines); (Inset) Single realizations in the 
Bloch sphere; measurement time Fr = 0.1, feedback param- 
eter: Qt = -0.45. Parameters: e = 0, U/T = 1, F'/F = 9, 
Tc/F = 1. 



arc simultaneously coupled back into Tc. This yields a 
frequency shift linear in Ot ^ ^lt = Qrt [compare 
Fig.mjb)] of the A-peak which is clearly revealed in S{ll!) 
[Fig. Efa)]. Moreover, at low frequencies the noise ex- 
hibits a pronounced local minimum w.r.t. Qt at small 
negative values (not shown) which turns out to corre- 
spond to the stabilization of pure qubit states. 
Generation of pure states - The simulation of the quan- 
tum trajectory in Fig.[3Ub) indeed confirms that by start- 
ing in a complete statistical mixture, the qubit will finally 
approach the pure state ^^jd*) — * 1^)) {S^ = —1) after 
a few electron jumps through the SET. In the following 
we discuss in more detail the stabilization of such pure 
qubit states. The steady state yields 

2 ' ^ ATciATc + sG) + c-f+{G + 2??^) ' ^ ' 

with s = sin(20T), c = 1 — cos(28t), and 
G= ^iVr + Vr)^. The 5^-component is solely gov- 
erned by electron-phonon relaxation [ll| (we omit the 
5^-component for brevity). We note that the SET oc- 
cupation is not affected by the feedback due to the high 
SET bias hmit. The steady-state purity [P = Tr(p2) = 
^^(S*")^], i.e., the length of the Bloch vector for smaU 
electron-phonon coupling 77^ <C F reads 



P={j-fc 



cG^ + 8Tc{4Tc + sG) 
c-f+G + 4Tc{4Tc + sGy 



(6) 



First, for zero feedback {Qt = 0, c = 0) the purity 
vanishes as expected. Second, above a critical cou- 
pling F'/F > 1 -I- 8Tc/T following from the condition 
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FIG. 4: (color online), (a) Location of purity P maxima in 
feedback parameter plane {Qt, T'/F), Eq. ((7|). Above a crit- 
ical coupling r'/F > 1 + 9>Tc/T (bifurcation) the qubit state 
is pure with P = 1, otherwise P < 1. At the bifurcation the 
qubit state is maximally delocalized with 5''' = —1 {tp = ti/2) 
and above / with Q < (p < it, tp ^ it /2 holds. The 
bifurcation becomes blurred and the complete delocalization 
is diminished for U/F = 0.5 (dashed line) and/or rjr/T = 
0.1, /3Tc — 1 (dotted line), (b)-(d) Normalized steady-state 
Current I/Io with Jo = (red, solid curves), Fano fac- 

tor F = 5'(0)/(2J) (blue, dashed curves), purity P = Tr(p^) 
(black, dotted curves) vs. Ot- Parameters: e — U — 0, 
n/F = 1, F'/F = 4 (b), 9 (c), 15 (d). 



B — \/(7 )^ — 167)? > 0, two maxima of P exist: 



e 



G[(7+)^ i:> + 12r,^] - 2-f+D ± 8T^B 
G{rT' + 16T2) 



(7) 



with D = + where the stationary qubit 

state has purity P = 1. Below the bifurcation 
{B^ < 0) the state is partially mixed (P < 1) and 
cos (St) = -G{G^ + 16T^)-^/^ holds. Fig. HJa) shows 
this bifurcation behavior of the maximal purity for var- 
ious Tc- The bifurcation point (B ~ 0) corresponds to 
the completely delocalized state with = —1 {(p = 7r/2) 
The sign of is determined by 7~, compare Eq. ([5|). In 
the upper(lower) branch the pure qubit state becomes 
increasingly localized with > (S"^ < 0). The ef- 
fect of finite U and the phononic environment onto the 
steady-state purity is essentially a blurring of the bifurca- 
tion in Fig.|4lja) (dashed and dotted curves). A complete 
delocalization {S"^ — —1) as for the noninteracting case 
hence cannot be obtained here anymore. Nevertheless, 
asymptotic purification with increasing localization can 
be realized in the limit of large SET-qubit coupling. 
Steady-state SET transport. - Without feedback we ob- 
tain the steady-state current Iq/b = (F -|- F')/8 and the 

Fano factor p = ^ = i + 2^^f±^^id^^^. The SET 



current is not affected by the phonons, whereas the noise 
carries information about the interaction of the qubit 
with the SET detector 75 and the phononic environment 
ijp. The Fano factor is one-half for vanishing SET-qubit 
coupling (F = F') and diverges for large SET-qubit cou- 
pling (F'/F ^ 1) which results from the bistable cur- 



rent behavior (quantum Zeno effect). It also diverges for 
Tc — ?> since the two subspaces decouple in this limit. 
Figs.|31Jb)-(d) present the Br-dependencies of the steady- 
state current and the Fano factor at various detector- 
qubit couplings F'/F. The current, in principle, follows 
the electron localization in the qubit. At the critical cou- 
pling [Fig.|4l^c)] when P = 1 the qubit state is maximally 
delocalized so that the current approaches its no-feedback 
value Jo and the Fano factor displays a local maximum. 
For larger couplings [Fig. W[d)] the purity maximum in 
the upper branch is accompanied with a current maxi- 
mum and a Fano factor minimum. Hence, the degree of 
the qubit state purity can be very accurately monitored 
by varying the feedback parameter Qt- 
Conclusions. - We have presented a simple way to sta- 
bilize exact pure states and generate charge oscillations 
in an individual charge qubit by applying an electronic 
feedback loop. Even for moderate detector back-action 
and coupling to a dissipative environment (e.g. phonons) 
an asymptotic purification can be achieved. We note that 
in contrast to ground state relaxation just by cooling and 
strong electron-phonon coupling our method also works 
at high temperatures and, in principle, arbitrary pure 
states can be produced. The detector current and its 
noise spectrum are shown to be sensitive monitors for 
the effect of feedback. 
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Equation numbers in the main text are referred to by an asterisk. 



QUBIT DYNAMICS WITHOUT FEEDBACK 



We rewrite Eqns. (2*) as vector equation x = C • x + v with 



c = 



— T 

-1l % 



— T 

1r 



no Ao+Vl Br 
ni Bl Ai + Vr 

X EE (n,So,Si)^, VEEE (7+,7+e^+7^e^,-7^e;j)^. 

I -1 



For its solution we use the Laplace transform x(z) = [zl — C] -(xq + v/z 

V = Irn Q2;0 qvO czO qxQ qyO qzO\T 

For symmetric coupling (r/2 = F^ = Tr, V /2 = F'^ = F'^) and e = [/ = we find 



(1) 



with the initial state 



Ai{z,no) + zA2{z) + (7+ + 2noz)A3{z) 
2z[A,izA) + {j+ + z)A,{z)] 



sy{z) = 2 



27^5*^0 + zSy'^ 



g{zy ' - - 8T^ + zg{z) 

ATcSy° + g(z)g-'° 
8T2 + zgiz) ' 



(2) 



with Ai{z,no) = Vrrirr + 2^02(7+ + z)] - (7^)^(7+ + 2z), ^2(2) = ^/^^7+ - 7-[.g(z)S'^- - ATcS'y-], 
Asiz) EE 8T2 + (7+ + z)(7+ + 2z), Sy^'- 



.y/zO 



3^^^'-°, and g{z) ee 7+ - VFF + 2z. 



The inverse Laplace transform for e = t/ =0 yields 



7s + 277^coth (/?T,) 



lint 



5*^/^° cosh (cjgi) =F — sinh (wgt) 



(3) 



with ujQ = i-\/7|) — (STc)^ (for 7^1 ^ 8Tc damped qubit oscillations occur.) and 7d = 75 + ^t]^ where 75 = 

■i(\/r — \fV'^ is the decoherence rate due to the SET-qubit coupling which was also found for a qubit coupled to a 
QPC in Ref. [ij where the respective rates are given by the QPC transmissions T and T' . For ?7 > the expressions 
become rather cumbersome such that the decay can be non-exponential (see e.g. Ref. and discussions therein). 
The steady state can be readily obtained by lim2^o+ zx(z) and yields (also for e, U 7^0) 



-^1 g _ / 27?., 

2' V7s + 2ry,coth(/3rc) 



0,0 



(4) 



with r\j. = girTcC^'^'^"/^''- being the phonon relaxation rate where g is the dimensionless electron-phonon coupling 
strength and Wc is a cut-off frequency 0]. For zero phonon temperature (/? — )■ 00) the .t— component becomes 
y/{l + y) {y = 2r]r/js) also discussed in Ref. the qubit relaxes into a pure state for y 3> 1 (qubit ground state). 
For large temperatures (/3 — 0) the qubit will behave as statistical mixture. 
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FEEDBACK-MODIFIED EQUATIONS OF MOTION 



The feedback scheme modifies Eqns. (2*) such that Ba — > B[, Uq — >• Hq, rii — > n(, and v — > v-^. For QaT = we 
find tIq = riQ, n( = rii, = v, and 



COs(eae)yr\jT Sm{&ae)V^a 

6f = I -sin(eae)yr^ cos(eae)yrj^ o | , 





(5) 



7^ 



which represents a rotation of the Bloch vector around the 2— axis with angle &ae- Respectively for Qae = the 
matrix B[ rotates the Bloch vector around the x— axis with angle QaT- 



Bl = 



v/Tji: 

COs(2eaT)v/r^jT sin(2eaT)7a+ 

-sin(2eaT)Vrj;^ cos(2e,T)7^ 
-7^(^0,sin(2eflT),cos(2ei^T) + ^) , 

7^ 0, sin {2Qlt), cos {2&lt) + ^ 
V 1l 

7l -7l [sin {2QLT)ey + cos {2QLT)ez]^ 



(6) 



For symmetric coupling and e ~ U ~ the a;— components decouple, and we can truncate the dynamical system to 
5 equations: 



V \{1 + cl) 



-\-np 2T, 

-2Tr. -4- 



/TT' 7^ 



SL -^CL 







21 \ 
2 ' 



/TT' 7^ 

•Jtt' 7''" 

^^Sfl ^CR 



2T, 



(7) 



with Sa — sin {2QaT) and Ca = cos (28aT) 
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